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Many electromagnetic problems involving more than one dielectric medium are not 
susceptible of an exact solution, when the appropriate boundary conditions are considered . 
The purpose of the present paper is to formulate a new boundary condition, which is capa- 
ble of leading to mathematically tractable problems, with limited sacrifices in accuracy. 



1. Introduction 

The theoretical treatments of mixed-path propa- 
gation [Wait, 1962a] problems usually depend upon 
the introduction of a Leontovich, or impedance type 
boundary condition which specifies the ratio of 

tangential E and H to characterize ground conduc- 
tivity. Otherwise, if one attempts to match tangen- 

-> -^ 

tial E and H across an interface, as required by the 
rigorous theory, then intractable boundary value 
problems usually arise. As a result, the use of linear 
boundary conditions to model the physics at an inter- 
lace is virtually required, if numerical results are 
desired. However, the Leontovich boundary condi- 
tion is known to represent the phenomena well only 
if the ground conductivity is high, in which case the 
surface impedance is reasonably independent of the 
angle of incidence. That is, if the ground has low 
losses, or is a fairly good dielectric, then there is 
significant penetration into this medium, and higher 
order boundary conditions are required if the physics 
at the interface is to be modeled accurately. It is the 
purpose of this note to introduce a more accu- 
rate version of the Leontovich boundary condition 
for use in propagation problems involving dielectric 
interfaces. 



2. Construction of the Boundary Condition 

In our search for a boundary condition that repro- 
duces the phenomenology of transition conditions at 
dielectric interfaces we have been motivated by the 
form of solutions to some radiation problems. The 
standard two-media problem of a radiating line source 
above a dielectric interface yields a secondary field or 
diffraction contribution, represented by an integral 
of the form 
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in addition to the free space Green's function. In 
the integral 7i(0) is the appropriate Fresnel re- 
flection coefficient for the dielectric interface, and 
I he contour (7 is the familiar path defining the Hankel 
function. The diffracted field has the interpretation 
of a summation of plane waves, traveling in all 
directions, real and imaginary, which have as a 
weight factor the Fresnel reflection coefficient extended 
to the complex 0-plane. This suggests that the 
scattered field in the dielectric half space contain- 
ing the source is characterized by the interface's 
reflection properties. 
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Figure 1 . The reflected fields for both the two-media problem 
and its reformulated counterpart have the same functional 
form and differ only iin the value of the reflection coefficient R. 

When medium k 2 is almost conducting it is known 
that an impedance type boundary condition of the 
form [Grunberg, 1943; Leontovich, 1944] 



du 
dy 



+\u=Q, 



(1) 



represents the phenomena fairly well. Therefore, 
for a more involved situation we seek to construct a 
higher order boundary condition with additional co- 
efficients so as to match the Fresnel reflection coeffi- 
cient. We start by noting that the boundary condition 



ih ^ +Au+ M 5? u= 



0, 
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(2) 
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imposed upon the wave equation 

v 2 u+kiu=o, y>o,h=-> (3) 

c 
implies plane wave solutions of the form (fig. 1) 

u=u iQC +R A (d)u te f, 
where 

„, ^ffciCxcos 6— V sin 0) 

^inc — & ? 

„, ^z/fciGzcos 0+2/ sin 0) 

W'ref — C j 

provided that the reflection coefficient R A {6) is chosen 
as 

cos0-(A-.£sin 2 0) 



Ra(P) = 



cos d+(A-B sin 2 6) 



(4) 



An approximation procedure is evident: find A 
and B such that R A {6) is a suitable approximation 
to the Fresnel reflection coefficient appropriate for 
the polarization of the excitation field. We em- 
phasize that the reflection coefficient in (4) is an 
approximating version of the Fresnel reflection 
coefficient by using the subscript A. For transverse 
electric excitation u(x, y) represents E z , E x =E y =0; 
and the correct Fresnel coefficient is [Stratton, 
1941], 



cos v— n 



i?TE (0) =" 



cos 6-\-n 



1- 4sin 2 



(5) 



The geometry is that of figure 1, and n=k 2 /k 1 , is the 
index of refraction where k 2 is the propagation 
constant in the lower half space. For transverse 
magnetic excitation u(x, y)=H 2) H x =H y =0] the 
reflection coefficient is 



Rtm(0)~- 



cos 6 

n 



cos 6-\ — 

n 



1 ^ snr 



1 k sin 2 1 



(6) 



For convenience only, we have assumed the magnetic 
permeabilities ^ and \x 2 of the two media to be 
identical. A comparison of 2? T e(0) and Rtm(Q) with 
(4) shows that we can approximate them for all 
real angles of incidence by picking A and B properly. 
To make this selection, we choose 
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Graph No. 1. 



The absolute percentage errors and the absolute errors that arise between the Fresnel reflection coefficients (5) and (6) , 

and their approximations (4) • 

The curves are plotted as a function for for the index of refraction n=l.Q. In this case B(d e ) has been chosen to make the match exact for 
e = O°, or normal incidence. 
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so that the approximation (4) will agree exactly 
with either (5) or (6) for normal incidence (6=0). 
For glancing incidence (0=±7r/2), the expression 
(4) will agree with (5) or (6) regardless of the choice 
of A or B since the forms of (4), (5), and (6) are 
such that all are equal to —1 for d= ±tt/2. 

Having motivated a choice of the constant A, we 
shall devote the next few subsections to selection of 
the remaining constant B, which can be chosen in 
many useful ways. 

2.1. General Angular Matching 

The constant B in (4) above may be used to make 
R A (0) exact for some additional angle e , where 
O<0 c <7r/2. For any choice of 8 ej the -Btm.te^'s 
are determined by matching 



1- 



-^S*™ sin 2 6 e with [~l-- 2 sin 2 d e T 

^TE, TM L W J 



where n = k 2 /ki is the index of refraction; inverting 
for B TETM we find 



£t 



!&)= 



Ate, tm { 1- I 1-^2 sin2 6 e\ J 

sin 2 B e 



(8) 



Observe that choosing 6 e is equivalent to specifying 

^TE. TM- 



2.2. Normal Matching 

If 7&>1, we can expand the radical in (8) by the 
binomial theorem, and use the A's defined by (7) 
to obtain 



*™=»[£r-i}F ain "'' + ---] 
^[^-|^ sin2 ^ + ---] 



Btm — 



(9) 



These leading terms of B TE and B TM are independent 
of 6 e the fitting angle; for analytical purposes it may 
be convenient to choose the B's as 



#i 



1_ 

~'2n 



B ™~2ri? 



(10) 



This selection of B will be referred to as normal 
matching since it corresponds to improving the agree- 
ment of the reflection coefficients in the neighbor- 
hood of 0=0, or normal incidence. The accuracy 
of this method is illustrated in graph la, where we 
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(a) 



Graph No. 2. 



(b) 



The absolute percentage errors and the absolute errors that arise between the Fresnel reflection coefficients (5), (6) 
and their approximations (4) plotted as a function of 6 for the index of refraction n = 1.6. 
In this case Ji(0,) 1ms been chosen to make the match exact for Brewster's angle 0e=0 6 =arc tan n, or for this case 58°. 
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plot the absolute percentage error ' PE(0) 



PE(0) = 



Ra(8) — ^te,tm(*) 



i? 7 



.w 



xioo 



(11) 



between the reflection coefficients (5) and (6) and 
their approximations (4) with normal match coeffi- 
cients for 7i=1.6. In the neighborhood of the 
Brewster's angle B =arc tan n, for which R T m(O b ) = 0, 
the percentage error, but not the absolute error, for the 
TM approximation is necessarily unbounded. In 
graph lb, we inspect the absolute error AE(0) as a 
function of 6 for the same fitting 



AE(d) = \R A (d)-R TEtTM (d)\ 



(12) 



and see that the absolute errors for both the TE and 
TM approximations have about the same values. 

2.3. Brewster's Angle Matching 

The unbounded percentage errors for the TM 
case in the vicinity of Brewster's angle can be 
eliminated by choosing B TM (d e ) to yield a perfect 
match at Brewster's angle 6 e —B B) and (8) becomes 



^te.tm(^)- ^te,tm^2~ i 1— 1— 2_j_ | r 



(13) 



Graphs 2a and 2b illustrate the percentage and 
absolute errors for this choice of B(6 B ) for n= 1.6. 
Notice the appreciable reduction in error as compared 
to the situation described by graph 1. 

2.4. Chebyshev Matchings 

Another choice for B is possible, at least for the 
TE approximation; this selection (fig. 2a) is a value 
of 5 T e so chosen as to minimize the maximum absolute 
percentage error for the TE reflection coefficient. 
The magnitude of the error changes sign when the 
match is exact, i.e., at 6=d e . Suppose we vary 
6 e until we obtain a value 6 e =6 T such that the maxi- 
mum negative percentage error was the same mag- 
nitude as the maximum positive error. This balance 
will minimize the maximum absolute percentage error 
over the interval —t/2<6<tt/2, by the Chebyshev 
criterion. The mathematical analysis to find the 
optimum Chebyshev matching would be interesting. 
Here we have done it by a trial and error process; 



1 We inspect errors rather than compare actual and approximate curves since 
the two graphs would appear to coincide if superimposed. 





Figure 2. (a) A Chebeyshev fitting is obtained when the fitting 
angle 6 e is so chosen that the maximum positive error, -fE max , 
equals the maximum negative error, — E max . (b) A two-point 
Chebeyshev fitting can be obtained by sacrificing the demand 
for a perfect match at = 0. 
In case (b) the maximum error emax would be less than the Emax of part (a). 

the values of 9 T required to afford this match for 
the TE case are plotted in graph 3. The maximum. 
associated errors for any real 6 in the approximate 
reflection coefficient (4) are plotted in graph 4 as a 
function of n. It is not possible to repeat this 
procedure for the TM case due to the unbounded 
nature of the percentage error in the vicinity of 
Brewster's angle. However, a Chebyshev matching 
is feasible for both the TE and TM cases by minimiz- 
ing the absolute error, but this was not analyzed. 
Better Chebyshev fittings (fig. 2b) might arise by 
relaxing the demand of a perfect fit at — 0; this per- 
mits choosing both A TE and B TE to improve the 
match. This procedure would yield a perfect fit 
for six values of d over the interval — ir/2 <&<-{- ir/2. 
This particular matching was not carried out 
because of the labor involved in doing it empirically, 
and the excellent fit of the initial approximations. 



2.5. Improved Angular Matchings 

For other problems a useful procedure is to im- 
prove the approximation in the neighborhood of 
a specified angle of incidence. An example of a 
situation that would motivate this procedure is 
a transmitter T and distant receiver R both located 
near the earth's surface. The geometry is such 
that much of the energy collected by the receiver is 
associated with the arrival cf rays in the neighbor- 
hood of grazing incidence. For this problem one 
might choose a value of B to match the derivatives 2 
of the exact and the approximate reflections coeffi- 
cient at grazing incidence; a simple calculation 
yields the result 

5 TE(TM (7r/2)=^ TE)TM i 1 "^-^] }' (14) 



2 The values of these reflection coefficients themselves already agree for 0=- 
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The corresponding values of B(8t) are given by (8). 



2.6. General Recommendations 

Unless we wisli to emphasize accuracy for a 
specified direction, it is reasonable to use the bound- 
ary condition (4) with Brewster's angle coefficients 
for problems of transverse magnetic excitation. For 
a problem of transverse electric excitation we would 
use the values of B T in curve (3) in formula (8) to 
yield a Chebyshev fit. The maximum percentage 
errors that arise for real angles of incidence for these 
two approximations are plotted for reference in 
graph 4 as a function of n, the index of refraction. 



3. Generalizations 

In addition to the above-mentioned methods, other 
refinements are possible. For example, any desired 
degree of accuracy can be obtained by adding addi- 
tional terms to the boundary condition: an inclusion 
of additional terms proportional to 

An the boundary condition (2) would have generated 
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Raw- 



cos 6-{A+B sin 2 fl+(?sin 4 0+D sin 6 6+ . . .) 
= cos d+(A+B sin 2 0+(?sin 4 d+D sin 6 6+ . . .) 



as the approximate reflection coefficient. The coeffi- 
cients of these higher order terms would then be 
available to match either as many terms in a series 
expansion of the radical 

fl-ism^T =A+B sin 2 S+Csm'd+D sin" 6+ . . . 

as needed for extreme accuracy, or for use in more 
refined Chebyshev approximations. 

From the preceding analysis, it is clear that the 
modeling of an interface by a linear boundary condi- 
tion need not be restricted to dielectrics. The pro- 
cedure can be used whenever the reflection coefficient 
from some arbitrary medium is given. For example, 
the reflection coefficient from an exponentially strati- 
fied medium is known [Wait, 1962b], and the present 
procedure can be repeated for that case. Such an 
analysis would be very interesting, and might be 
used, for example, to discuss ionospheric propagation 
from a day to a night region. 
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Figure 3. The geometry of a test problem which serves to compare 
the use of the new boundary condition (2) with the use of the 

-» -> 

continuity of tangential E and H at a dielectric interface in 
a problem involving a horizontal discontinuity. 



4. Utility of the Boundary Condition in a 
Problem Involving a Horizontal Discontinuity 

We have lent plausibility to a procedure which 
characterizes transition conditions at a dielectric 
interface by a linear boundary condition. This 
procedure is valid if the interface separates two half 
spaces, but what is its utility in a more complex 
geometry, say one for which there is a discontinuity 
in the horizontal direction? To answer this question, 
we have investigated a rigorously solvable test 
problem of the Wiener-Hopf type (cf , fig. 3) . A plane 
wave is launched along a semi-infinite metallic screen 
located at the interface of two dielectrics. The 
termination of the screen produces a diffracted field 
that depends upon the physical parameters in a very 
complex fashion. It will be very significant if the 
new boundary condition implies results comparable 
to that obtained by a rigorous matching of tangential 

E and H. 

We obtain a reformulated version by removing the 
lower medium k 2 and replacing it by the new bound- 
ary condition (2) at y=0, z>0. The coefficients in 
this boundary condition are chosen to characterize 
the former transition conditions at the k x —k 2 inter- 
face. We have indicated that there exists some 
flexibility in the choice of the parameters A and B; 
consequently, we shall leave them unspecified to see 
how various choices affect the diffracted field. 

Let us specialize the discussion for the TM-case 
for which H x =H y =0, since this polarization allows a 
plane wave to propagate along a perfect conductor at 
grazing incidence. The importance of this fact is 
that if the plane wave travels in the positive ^-direc- 
tion, we can compare the diffracted fields for £-^ + oo ? 
y=0 in both problems without the necessity of 
performing the explicit Wiener-Hopf decomposition. 
It would be of value to compare the fields in other 
directions, but that would involve a lengthy numeri- 
cal program. 



Let u(x, y) =H Z , then a standard approach [Kane 
and Karp, 1960], using dual integral equations, and 
Van der Waerden's saddle point analysis leads to the 
following asymptotic developments for x-^> + oo . For 
the exact problem, using matching conditions at the 
interface, we find for the far-field 



■&-*■•>— (k)'Gri) 



1/2 6 z(fc 1 x- 7 r/4) 



M 1 



-o(x- 3/2 ) 

(15) 



and likewise a corresponding result for the reformu- 
lated counterpart u A (x,y) using the linear boundary 
condition (2) along the positive £-axis 



\imu A (x,0)-- 



(2& ) 1/2 0*^1*- *-/4) 



A-B hiirx) 1 



0(x~ 3/2 ). (16) 



In both the rigorous (15) and the reformulated 
(16) asymptotic developments, the leading term in 
the diffracted field has canceled the incident plane 
wave. This is of course a consequence of the fact 
that the reflection coefficient for grazing incidence 
is —1 for either problem. It is more important to 
notice that for both problems, if x^>>1, the junc- 
tional jorm of the leading terms (15) and (16) agree, 
and dvffer only by a numerical factor 



k\ 



1 



1 



Jd (kt-k^ h(A TM - 
or, if we choose A TM =l/n, 

n , 1 
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kill 
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\n 



-B, 



The percentage error PE between these two coeffi- 
cients is 



PE= 



1— nB T y 



XI 00, 



(17) 



and, by (9), B TM is of the order 0(1 /n 3 ) for fitting at 
any real angle, so that this error decreases with in- 
creasing n. 

In line with the remarks in section 2.5, we observe 
that if in the far-field, the preferred direction had 
been chosen as grazing incidence, we would have 
used (14) to select a B TM to match the derivative of 
the reflection coefficient at grazing: incidence, i.e., 



#1 



-o-# 



(18) 



For this choice there would have been no error at all 
in (15) as compared with (16) for any n. 

We emphasize that the functional form of the lead- 
ing terms in the diffracted field are the same for 
y= 0, #>>1, but with differing coefficients for var- 
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ious choices of B TM . For any selection of this parameter 
we can compute the percentage error (17) between 
the leading coefficients of (15) and (16). In graph 
5, we plot three curves: (A) represents the "normal 
impedance" or Leontovich type approximation 

-tj-- r- — \-A TM u=0, obtained by setting B TM = 0; (B) 

shows the error obtained by using the normal match 
values for B TM given by (9). The smallness of the 
error in this case is significant because we are looking 
at the diffracted field in a direction perpendicular to 
that for which the coefficient B TM was selected. 
The curve (C) illustrates the error if B TM is chosen 
as indicated by (13) which is a matching at Brewster's 
angle. We obtain the expected improvement over 
case (B) since Brewster's angle is closer to grazing 
incidence than normal incidence. Notice that the 
big improvement is from (A) to (B), a correction 
which arises from the additional term in the boundary 
condition. 

5. Concluding Remarks 

5.1. Reciprocity 

A feature of the new boundary condition (2) is 
that it involves a second derivative. This somewhat 
surprising detail has been inserted to guarantee 
reciprocity. A first derivative could have been 
used, but it would not lead to reciprocal phenomena. 
This point can be most easily established by an 
inspection of the associated reflection coefficient. 
If the boundary condition had been written with a 
single x-derivative 

1 bu , ,Bdu_ 
iki dy ki dx 

its reflection coefficient would have been 



RaW-- 



cosfl— (A— B sin 6) 
"cos 0+ (A— I? sin 6) 



which is not the same for a wave reflected from the 
left as compared to one reflected from the right since 
sin 6 changes sign with 6. 

5.2. Existence and Uniqueness 

That solutions of problems involving the new 
boundary condition do exist can readily be shown by 
construction [Kane and Karp, 1960; Jones, 1962]. 
In addition these references illustrate the use of the 
boundary condition, and also its ability to reproduce 
complex diffraction phenomena. In addition, the 
uniqueness has also been established [Kane, 1961]. 

5.3. Limitations 

The approximation procedure is limited by the 
restriction k 2 ^>k h or n>l. For, if k 2 <^k h there 
arises a real angle of incidence corresponding to total 
reflection at which the second radical in (5) and (6) 
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Graph No. 5. 

A comparison of the errors that arise for x»l,y = 0m the total 
field including diffraction effects of the test problem involving 
a horizontal disco ntinuity. 



becomes imaginary. Since (4) involves constant, 
coefficients, the boundary condition is unable to 
imitate a transition from a real reflection coefficient 
to a complex one for real angles of incidence. This 
is clear by inspection of the reflection coefficient (4) ; 
so long as both A and B are real, the reflection 
coefficient remains real for any real angle 6. 

Observe, from graph 5, that the approximation 
deteriorates for ?i~l. Physically n~l corresponds 
to an interface with a vanishing reflection coefficient, 
which is difficult to approximate on a percentage 
basis. 4 However, the absolute errors remain small 
and bounded. As a consequence the boundary 
condition (4) specialized for n=l may be useful 
as a model of a perfect absorber since the reflection 
coefficient it defines is so small. A curve of this 
special R A (0) is drawn in graph 6. We see that 
R(6)^0 until |0|=6O°;it is unavoidable that 7i(0)-> 

— 1 for |0|->k> this is so even for the exact reflection 

coefficient if n^\. 



* However this fact allows one to use other approximate methods [Karp and 
Sollfrey, 1950]. 



309 



LU 
O 

o 

z: 
o 

h- 
o 

QJ 































^MAX = .02 
































1 

i k 


du + A + _B_ iu 2 _ 














^ 1.0 
) = 0.6 


* 
















E 


414 


















1 













0° 10° 20° 30° 40° 50° 60° 70° 80° 90° 

ANGLE OF INCIDENCE 6 

Graph Xo. 6. 

A special boundary condition which defines an almost vanishing 
reflection coefficient for < |0| <60°, indicating that the boundary 
condition may be used as a model of a perfect absorber. It 
is unavoidable that 7?(0)— > — ./ for \9\=t/%; this is forced by 
the form of the reflection coefficient {2.8) . 
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